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Impact of nuclear effects on the determination of the nucleon axial mass 
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We analize the influence of nuclear effects on the determination of the nucleon axial mass from 
nuclear cross sections. Our work is based on a formalism widely applied to describe electron-nucleus 
scattering data in the impulse approximation regime. The results of numerical calculations show 
that correlation effects, not taken into account by the relativistic Fermi gas model, sizably affect the 
Q^-dependence of the cross section. However, their inclusion does not appear to explain the large 
values of the axial mass recently reported by the K2K and MiniBooNE collaborations. 
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Experimental searches of neutrino oscillations exploit 
neutrino-nucleus interactions to detect the beam parti- 
cles, whose properties are largely unknown. The use of 
nuclear targets as detectors, while allowing for a substan- 
tial increase of the event rate, entails non trivial prob- 
lems, since data analysis requires a quantitative under- 
standing of neutrino-nucleus interactions. In view of the 
present experimental accuracy, the treatment of nuclear 
effect is in fact regarded as one of the main sources of 
systematic uncertainty (see, e.g., Ref.[l|). 

The description of nuclear dynamics is particularly 
critical to analyses aimed at obtaining nucleon proper- 
ties from nuclear cross sections. 

Recently, the K2K and MiniBooNE Q collabora- 
tions have determined the nucleon axial mass Ma, i-C. 
the mass scale driving the Q^-dependence of the dipole 
parametrization of the nucleon axial form factor, from 
neutrino interactions with oxygen and carbon, respec- 
tively. The reported value, Ma ~ 1.2 GeV, turn out 
to be significantly larger than the one previously deter- 
mined from deuterium cross sections. Ma ~ 1.0 GeV [1]. 
The authors of Ref . 0] argue that the large Ma extracted 
from the data should be regarded as an "effective ax- 
ial mass" , embodying nuclear effects not included in the 
Relativistic Fermi Gas (RFG) model employed in their 
analysis. They also suggest that replacing the RFG with 
one of the more advanced nuclear models available in the 
literature [1, H, 0, Q may result in a value of Ma closer 
to that measured using deuterium. 

This letter is aimed at assessing the impact of the treat- 
ment of nuclear effects on the determination of the nu- 
cleon axial mass. Our work is based on the approach 
described in Refs.j^, [^, in which nucleon- nucleon corre- 
lations not included in the RFG model are consistently 
taken into account. 

Both K2K and MiniBooNE search for signatures of 
neutrino oscillations using Charged Current Quasi Elas- 
tic (CCQE) interactions 



the cross section at neutrino energy < 1.5 GeV [l^l- 

The differential cross section of process ([T]), in which a 
neutrino carrying four- momentum = [E^, \<i^) scatters 
off a nuclear target producing a muon of four-momentum 
/c^ = (£"^1 k^), while the target final state is undetected, 
can be written in Born approximation as 
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where Gf is the Fermi constant and Vud is the CKM 
matrix element coupling u and d quarks. The tensor 
Lap is fully specified by the lepton kinematical variables, 
while the definition of W'^^ involves the target initial and 
final states, as well as the nuclear weak current. 

For neutrino energies larger than ^0.5 GeV, can 
be obtained within the impulse approximation (lA), i.e. 
assuming that neutrino-nucleus scattering reduces to the 
incoherent sum of scattering processes involving individ- 
ual neutrons, whose momentum (p„) and removal energy 
{E) distribution is described by the nuclear spectral func- 
tion P(p„,i?) [ll|,[l3]. Neglecting final state interactions 
(FSI) between the struck nucleon and the spectator par- 
ticles, the nuclear tensor can be written in the form [^, Q 
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is the four-momentum transfer and 



where q = k 

p = (i?„, p„), with En = ("4 + |PnP)^^^, rUn being the 
neutron mass. 

The tensor W"^ describes the charged current weak 
interactions of a neutron of initial momentum p„ in free 
space. The effect of nuclear binding is accounted for by 
the replacement q ^ q = (qo, q), with 



<?o = 90 + Ma - Ea^i 



En 



(1) where Ma and Ea-i = \{Ma - to„ + E)^ 
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which are known to yield the dominant contribution to 



denote the target mass and the energy of the recoiling 
nucleus, respectively @. 
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The above procedure, originally proposed in the con- 
text of a study of electron induced nucleon knock out 
processes , accounts for the fact that a fraction of the 
energy transfer to the target goes into excitation energy 
of the spectator system. The energy Sqq = qa — qa is spent 
to put the struck particle on the mass shell, and the ele- 
mentary scattering process is described using free space 
kinematics with energy transfer qo- The physical intepre- 
tation of go emerges most clearly in the (|Pn|/77in) ~> 
limit, corresponding to qo ~ qo ^ E. 

The main effect of FSI on the differential inclusive 
cross section is a redistribution of the strength, result- 
ing from the coupling of the one particle-one hole final 
state to more complex n-particle n-hole configurations 
[l5l |. This leads to a sizable quenching of the cross sec- 
tion in the region of the quasi free peak, corresponding 
to qo ^ Q^/2mn, with = —q^, associated with the 
enhancement of the tails at both low and high qo- 

In this work FSI have been described following the ap- 
proach originally developed in Ref . . in which the in- 
clusive cross section is written in the convolution form 
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where {da-Q / dil ^dE'^) is the cross section in the absence 
of FSI , obtained from Eqs.((l]) and ([3]). The folding func- 
tion, embodying FSI effects, is trivially related to the 
spectral function of particle states [l^. It can be com- 
puted using the formalism of nuclear many-body theory 
and the eikonal approximation, i.e. assuming that: i) the 
outgoing proton moves along a straight trajectory with 
constant speed, and ii) the spectator nucleons act as a 
collection of fixed scattering centers [1, ■ 

The results of Ref. 0] show that inclusion of FSI effects 
is needed to reproduce the measured cross sections of 
the process e O ^ e' + X at electron beam energy 
- 1 GeV. 

In order to assess the validity of the analysis of Refs.Q, 
0] and the role of nuclear effects not taken into account l3y 
the RFG model, we have first studied the dependence on 
Ma of the Q^-distribution of the process - 
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p+ (A — l) at beam energy Ei, —1.2 GeV, corresponding 
to the peak in the energy spectrum of the neutrinos used 
by the K2K collaboration Q]. 

In Fig. [1] we compare the RFG results to those ob- 
tained using Eqs.([l]) and ([3]) and the spectral function of 
Ref. [Ill, with Ma = 1.0 and 1.2. 

It clearly appears that, while increasing the axial mass 
leads to an enhancement of the cross section, the in- 
clusion of correlation effects through the use of a real- 
istic spectral function produces a sizable quenching^The 
dashed line, corresponding to the model of RefsTjl, 
turns out to be below the dot-dash line, corresponding 
to the RFG model, over the whole range. Note that 
FSI effetcs, not taken into account in the calculations, 
lead to a further suppression of the curves labelled SF. 

Based on the results of Fig. [2 we conclude that, as far 
as the distribution at fixed neutrino energy is con- 
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FIG. 1: (Color online) Q'^-dependence of the cross section 



of the process 



O 



^ + p + X , for neutrino energy 



Ey = 1.2 GeV. The dashed and solid lines, labeled SF, have 
been obtained using the approach of Refs.[^, 0| with Ma = 
1.0 and 1.2 GeV, respectively. The corresponding results of 
the RFG model, with Fermi momentum pf = 225 MeV and 
removal energy e — 27 MeV, are represented by the dotted 
and dot-dash lines, labeled FG. 



cerned, a larger value of the axial mass cannot be ex- 
plained by replacing the RFG with the more advanced 
model of nuclear dynamics discussed in Refs.0, Q. 

It is important to realize, however, that using a realis- 
tic momentum and removal energy distribution may also 
significanltly affect the determination of E^. 

From the requirement that the elementary scattering 
process, i^^ + n ^ /i + X, be elastic, i.e. that 
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where is the proton mass and the four momentum 
of the struck nucleon is given by p„ = (i?„,p„), with 
En — Ma — Ea-1, it follows that 
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where 9^ is the muon angle relative to the neutrino beam 
and cos6'„ (k^ • p„)/(|k^||p„|). 

Setting |p„| =0 and fixing the neutron removal energy 
to a constant value e, i.e. setting E = e, implying in turn 
En ~ nin — e, Eq.([7|) reduces to 
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with Am^ = ml 



In the analysis of both K2K and 



MiniBooNE data, the energy of the incoming neutrino 
has been reconstructed using the above equation (com- 
pare to Eq.(5) of Ref.0 and Eq.(3) of Ref.[|), with e = 
27 and 34 MeV for oxygen Q and carbon 0], respectively. 
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Equationl!?!) clearly shows that, in general, the knowl- 
edge of and 0^, the observables measured by K2K 
and MiniBooNE, does not uniquely determine the neu- 
trino energy. For any given E^^ and 0^, E^, depends on 
both magnitude and direction of the neutron momentum 
Pn, as well as on its removal energy E, entering the defin- 
tion of En- 

The distribution of neutrino energy can be obtained 
from Eq.([2l) using values of |p„| and E sampled from 
the probabihty distribution |p„pP(p„,£') and assuming 
that the polar and azimuthal angles specifying the direc- 
tion of the neutron momentum be uniformly distributed. 
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FIG. 2: (Color online) Upper panel: Neutrino energy distri- 
bution at Efj, = 600 MeV and 6^ = 60°, reconstructed from 
Eq.l!?]) using 2 xlO^ pairs of (|p„|,i7) values sampled from 
the probability distributions associated with the oxygen spec- 
tral function of Ref.[l2| (SF) and the Fermi gas model, with 
Fermi momentum pf = 225 MeV and removal energy e = 27 
MeV(FG). The arrow points to the value of Ei, obtained from 
Eq. ((si . Lower panel: Same as the upper panel, but for iJ^ — 
1 GeV and 9^ = 35°. 

The spectral functions of nuclei ranging from car- 
bon to gold, computed in Ref . [l^ . exhibit high momen- 
tum and high removal energy tails, extending well above 
|p„| > 500 MeV and E > 100 MeV. A direct measure- 
ment of the carbon spectral function from the (e, e'p) 
cross section at missing momentum and energy up to 
^ 800 MeV and ^ 200 MeV, respectively, has been re- 
cently carried out at Jefferson Lab [l3| ■ The preliminary 
results of data analysis appear to be consistent with the 



theoretical predictions of Ref. [ij] . On the other hand, in 
the RFG model the typical Fermi momentum, pp, and 
average removal energy (e of Eq.®) are ~ 200 MeV and 
~ 30 MeV, respectively. 

To gauge the effect of the high momentum and high 
removal energy tails of P{pn, E), we have computed the 
neutrino energy distribution, F{E^), using 2 xlO^ pairs 
of (|p„|,i?) values drawn from the probability distribu- 
tions associated with the oxygen spectral functions of 
both Ref.[l^ and the RFG model, with Fermi momen- 
tum = 225 MeV and removal energy e = 27 MeV. The 
results corresponding to E'^ = 600 MeV and 0^ = 35°, 
and Efj^ = 1 GeV and 6*^ = 35°, are displayed in the 
upper and lower panels of Fig. [21 respectively. 

It appears that the distributions predicted by the RFG 
model arc more sharply peaked at the neutrino energy 
given by Eq.®. On the other hand, the F{E^) obtained 
from the spectral function of Ref. [iJI are shifted towards 
higher energy by ^ 20 MeV, with respect to the RFG 
results, and exhibit a tail extending to very large values 
oiE,. 

The average values of Ey calculated using the RFG dis- 
tribution turn out to be {E^) = 1.28 GeV and 934 Mcv 
for the kinematics of the upper and lower panels, respec- 
tively, to be compared with 1.28 GcV and 931 Mev re- 
sulting from Eq.®. On the other hand, the distributions 
associated with the spectral functions of Ref.|T2l yield 
{E„) = 1.35 GeV and 1.01 Gev. 

As Q'^ grows with E^ according to 



= 2E,E^ ( 1 - cos 
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the results of Fig. [T] and [2] suggest that using neutron 
energies and momenta obtained from a realistic spectral 
function in the equation determining the kinematics of 
the CCQE process would lead to extract an even larger 
value of the axial mass. 

The energy and momentum distribution of the struck 
neutron also affects the weak interaction vertex, since 
the neutrino interacts with a hound moving neutron. As 
pointed out above, within the lA formalism binding is 
taken into account through a shift of the energy transfer. 
As a result, within the approach of Rcfs.[^, Q the ele- 
mentary neutrino-neutron scattering process takes place 
at Q^ = H^-ql>Q\ 

In order to assess the full impact of replacing the RFG 
model with the approach of Refs. [1, , we have computed 
the differential cross section of the process 4- A — ^ 
^ + p -\- {A — 1), as a bmction of the incoming neutrino 
energy E'y, for the muon kinematics of Fig. [21 

The solid lines of Fig. [3l show the results of the full 
calculation, carried out using the spectral function of 
Ref. [l3|, while the dashed lines have been obtained ne- 
glecting the effects of FSI and the dot-dash lines cor- 
respond to the RFG model. Note that, unlike the his- 
tograms of Fig. [21 the curves displayed in Fig. [31 have 
different normalizations. As in Fig. [21 the arrows point 
to the values of E^, given by Eq.®. 
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FIG. 3: (Color online) Upper panel: Differential cross section 
of the process u^^+A ^ ^+p+{A—l), at = 600 MeV and 
9fj, = 60° , as a function of the incoming neutrino energy. The 
solid line shows the results of the full calculation, carried out 
within the approach of Refs. [1,01, whereas the dashed line 
has been obtained neglecting the effects of FSI. The dot-dash 
line corresponds to the RFG model with Fermi momentum 
Pf = 225 MeV and removal energy e — 27 MeV. The arrow 
points to the value of obtained from Eq. ([8]). Lower panel: 
Same as the upper panel, but for — 1 GeV and = 35°. 



The differences between the results of the approach of 
Refs. 0, Q and those of the RFG model appear to be 



sizable. The overall shift towards high energies and the 
tails at large E^, present in the histograms of Fig. [21 
are still clearly visible and comparable in size, while the 
quenching with respect to the RFG model is larger than 
in Fig. [21 Comparison between the dashed and dot-dash 
lines indicate that the ~ 20 % difference at the peak of the 
distributions is mainly due to the replacement —i- Q^, 
while inclusion of FSI effects leads to a further reduction 
of about 8%. 

In conclusion, the results discussed in this paper in- 
dicate that nuclear effects not included in the RFG 
model significantly affect the Q^-dependence of CCQE 
neutrino-nucleus interactions. However, contrary to the 
expectation of the authors of Ref. Q , their inclusion does 
not help to reconcile the large values of Ma reported 
in Refs.0, Q with those extracted from deuterium data. 
Using the model of Rcfs.j^, Q in the data analysis would 
in fact lead to predict an even larger value of the axial 
mass. 

Other possible explanations of the disagreement be- 
tween the values of Ma obtained by different experi- 
ments, such as misidentification of CCQE events, should 
be carefully investigated using state-of-the-art models of 
nuclear structure and dynamics. 

While this paper was being drafted, the NOMAD 
collaboration released the results of the analysis of 
quasi-elastic muon neutrino and antineutrino scattering 
data, yielding an axial mass Ma = 1-05 ± O.Q2{stat) ± 
OM{syst) GeV [Hi. As the NOMAD experiment takes 
data at beam energies much larger than those used in 
both K2K and MiniBooNE, we have chechek the sensi- 
tivity of our results to the incoming neutrino energy. It 
turns out thet the main features of the Q^-distribution 
of the process i^f^ +-^^ 0—^fi+p+{A — l) do not change 
significantly as E^, increases from 1.2 GeV to 10 GeV. 

The authors are indebted to M. Sakuda, for drawing 
their attention to the subject of this paper. Useful discus- 
sions with L. Ludovici, M.H. Shaevitz and M.O. Wascko 
are also gratefully acknowledged. 
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